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ABSTRACT 

On 1998 November 14, Saturn and its rings occulted the star GSC 0622-00345. We observed atmospheric 
immersion with NSFCAM at the National Aeronautics and Space Administration's Infrared Telescope Fa- 
cility on Mauna Kea, Hawaii. Immersion occuiTed at 55. °5 S planetocentric latitude. A 2.3 //m, CH4-band 
filter suppressed reflected sunlight. Atmospheric emersion and ring data were not successfully obtained. We 
describe our observation, light-curve production, and timing techniques, including improvements in aperture 
positioning, removal of telluric scintillation effects, and timing. Many of these techniques are known within 
the occultation community, but have not been described in the reviewed literature. We present a light curve 
whose signal-to-noise ratio per scale height is 267, among the best ground-based signals yet achieved, despite 
a disadvantage of up to 8 mag in the stellar flux compared to prior work. 
Subject headings: atmospheric effects - infrared: planetary systems - occultations - planets and satellites: 

individual (Saturn) - techniques: image processing - techniques: photometric 

Online-only material: tar.gz 



1. INTRODUCTION 

Stellar occultations by planets offer among the highest spa- 
tial resolution of any astronomical observing technique, com- 
peting even with atmospheric descent probes for certain mea- 
surements. From Earth, spatial resolution is limited by the 
projected stellar beam size and Fresnel diffraction, the latter 
giving ^1 km resolution for solar system giant planets for vis- 
ible and near-IR wavelengths. This is sufficient to reveal fine 
structure in a planetary ring or propagating gravity (buoyancy) 
waves in an atmosphere. However, the occultation technique 
is demanding, and failing to optimize even a single parameter 
of an observation can significantly reduce the quality of the 
derived light curve. Occultation observers have thus devel- 
oped numerous specialized techniques to improve their data. 
Many of these have never been described in the literature. 

On 1998 November 14, Saturn and its rings occulted GSC 
0622-00345, as predicted by Bosh & McDonald (1992). We 
obtained a light curve for atmospheric immersion, based on 
infrared imaging observations at the NASA Infrared Tele- 
scope Facility (IRTF) on Mauna Kea, HI. Figure 1 shows our 
viewing circumstances and Table 1 gives parameters of the 
event. 

Occultation cadences are typically from a few to 10 Hz, 
and in the case of lunar occultations they can be much higher. 
These rates are sufficient for the temporal resolution of flux 
variations caused by scintillation (variable refractive focus- 
ing and defocusing of the stellar beam) in Earth's atmosphere. 
This is typically the limiting noise source for occultations by 
bright stars. During our observations, we monitored Rhea, 
which experienced approximately the same scintillation ef- 
fects as the star, just arcseconds away. We were thus able to 
use measurements of Rhea's (assumed constant) flux to com- 
pensate for scintillation. This, in combination with other tech- 

^^ Visiting Astronomer at the Infrared Telescope Facility, which is operated 
by the University of JJawaii under Cooperative Agreement no. NCC 5-538 
with the National Aeronautics and Space Administration, Science Mission 
Directorate, Planetary Astronomy Program. 



niques, produced a light curve of very high quality. 

Here we present a full accounting of the methods we em- 
ployed to observe the event and derive the light curve, in- 
cluding those previously undocumented. Section 2 describes 
observing techniques, data reduction methods, and the light 
curve. It also discusses an attempt to apply optimal photom- 
etry to light curve production and thoughts on improving in- 
strument interfaces to avoid failures such as our lost egress 
observation. Section 3 presents a method for independently 
verifying the system timing solution. Section 4 offers brief 
conclusions. We include the light curve and software for ana- 
lyzing occultation timing data as electronic supplements. Pa- 
per II (Harrington et al. 2010) uses the derived light curve to 
investigate the atmosphere of Saturn. 

2. OBSERVATIONS AND LIGHT CURVE 

The NSFCAM InSb aiTay camera (Shure et al. 1994) 
recorded the event in "movieburst" mode. This mode writes a 
time sequence of images to computer memory and then saves 
the set to disk after the observations finish. We recorded three 
subarrays derived from the full 256 x 256 pixel array. These 
boxes contained the star (28 x 32 pixel), Rhea (28 x 32 
pixel), and a piece of blank sky (16 x 16 pixel), placed as 
shown in Figure 1. The larger boxes were originally 32 x 
32 pixels, but a camera software eiTor overwrote four edge 
columns. No critical data were lost and the error is now fixed. 
Table 2 gives additional observation and light-curve parame- 
ters. 

The limit to spatial resolution in a stellar occultation is gen- 
erally the larger of the Fresnel scale and the projected size of 
the occulted star. For typical in situ spacecraft occultations, 
both are very small compared to the sky-plane distance be- 
tween samples, but for Earth-based work either can dominate. 
In the case of 28 Sgr, the projected stellar diameter at the dis- 
tance of Saturn was ^20 km, much larger than the ^1-km 
Fresnel diffraction scale, resulting in significant smoothing of 
the ring and atmospheric occultation profiles (e.g., Harring- 
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TABLE 1 
OCCULTATION PARAMETERS 



Description 


Value 


Comment 


Star 


GSC 0622-00345 




R.A. 


Ih49m54"358 


J2000'' 


Decl. 


8°23'12"56 


J2000'' 


V magnitude 


11.002 


a 


K magnitude 


8.395 


b 


Spectral Type 


K5 




Proj. diam. at Saturn 


~0.7 km 




Half-light latitude 


55P5S 


Planetocentric 


Fresnel scale 


1.2 km 


VAd/2 


Sky-plane velocity 


17.940 km s' 


1 


Perpendicular velocity, vx 


12.802 km s^ 


vs. oblate limb 


Ring opening angle 


-14.° 6772 




Pole position angle 


l.°5134 


J20()0 


Saturn's geocent. dist., d 


1.25168x10^ km 





''Tycho-2 record 131066, H0g et al. (2000), apo, ~ 65 mas, V = 
VT -0.09Q{BT -VT), using catalog VT and BT magnitudes 
''2MASS 01495435-H0823123, Skrutskie et al. (2006) 



ton et al. 1993, Figure 4). To estimate the angular diameter of 
GSC 0622-00345, we first used the color index V-K = 2.607 
(Table 1) to infer a dwarf spectral type of K5 (Glass 1999). 
The stellar diameter for this spectral class (Lang 1992) yields 
an estimated angular diameter of 0. 12 mas, or 0.7 ± 0.2 km at 
Saturn's distance, where the error bar takes into account the 
possible effects of modest interstellar reddening. Even with 
a large uncertainty, the projected size of the occulted star is 
less than the Fresnel scale, which thus determines the intrin- 
sic spatial resolution of the occultation. 

The 2.3-/im filter is centered on a methane band, where Sat- 
urn's atmosphere strongly absorbs sunlight, so contrast be- 
tween the star and Saturn is very high. The rings are bright at 
this wavelength, but their spatial separation from the immer- 
sion latitude made it practical to subtract their scattered light 
(see below). 

At the half-light time (see below), seeing was 0."9 FWHM 
and the star was at 1.3 air masses and setting. The night was 
clear, but Rhea's flux varied by up to 8% in successive frames 
(2.3% rms) due to scintillation. Rhea was nearly as bright as 
the unocculted star (see Table 2) and had a point-spread func- 
tion (PSF) that very closely matched that of the star in each 
frame. We thus used Rhea as a standard for flux and position, 
and attempted to use it as a PSF standard. In contrast to the 
scintillation, sky emission was very steady and low, based on 
the separately recorded sky box. 

To obtain sub-pixel accuracy in image shifts, photometry 
aperture placement, and scattered-light subtraction, we sub- 
sampled by a factor of 10 using nearest-neighbor sampling (a 
factor of 5 gave essentially identical results). We performed 
all operations on the expanded images, except as noted be- 
low. These subpixels allowed a reasonable approximation to 
a circular photometry aperture that included partial pixels. 

After debiasing, flat-fielding, and interpolating three hot 
pixels that were not near the star or Rhea, we centroided Rhea 
in all images. The relative position of the star in two well- 
separated sets of 100 pre-occultation frames determined the 
Rhea-star offset and the offset's drift with time. This allowed 
accurate placement of the Saturn and ring background tem- 
plate (described below) and photometry apertures, even when 
the star was fully occulted. A third set of stellar centroids 
confirmed that the predicted stellar positions were accurate to 




+Star box 



Rhea box 




^H Sky box 

Fig. 1 . — Occultation geometry, subarray placement, and photometry aper- 
tures at the half-light time. The light line marks the path of GSC 0622-00345 
relative to Saturn, with crosses every 20 minutes. The heavier portion gives 
the span of the data set. The relative positions of the star, Rhea, and sky sub- 
arrays were constant. Telescope tracking kept the star in roughly the same 
pixel while Saturn and its moons drifted through. Tethys's motion relative 
to Saturn is evident by the position mismatch between the half-light data and 
the scattered-light template. We computed the template from post-occultation 
frames and subtracted it from the data to produce the star frame shown here. 



within 0."045 (rms), or 0.15 unsealed pixels. To estimate Sat- 
urn's position (for template placement, see below), we added 
the Rhea-Saturn relative sky-plane offset to Rhea's measured 
position in each frame. The Jet Propulsion Laboratory's Hori- 
zons On-Line Ephemeris System Web site^ generated the rel- 
evant positions on 2004 January 9 based on the DE-0406, LE- 
0406, and S ATI 36 ephemerides. Our derived Rhea-star drift 
rate closely matched the ephemeris prediction. 

The small amount of scattered light from Saturn and the 
rings complicated stellar photometry by requiring the subtrac- 
tion of a scattered-light template. We made a template by cal- 
culating the shifted, masked average (see below) of a set of 
post-immersion images, each divided by Rhea's flux in that 
frame and shifted to align the predicted position of Saturn in 
all frames. All shifts both in creation and application of the 
template were in units of integer subpixels, so no interpolation 
occurred. 

A shifted, masked average starts with two similar 3D 
"stacks," respectively containing the images and their bad 
pixel masks (in our case, all the masks are the same). In a 
mask, pixel values are one for good pixels and zero for bad 
pixels. The technique proceeds as follows: pad the stacks with 
zeros to prevent good data from shifting outside the bound- 
aries; multiply the image stack by the mask stack to set bad 
pixels to zero, shift the images to align them, and apply the 
same shifts to the masks; collapse each stack into a 2D image 
by summing along the image index dimension; at each pixel 
location, divide the image sum by the mask sum to calculate 
the shifted mask average. Each stacked image shifted differ- 
ently, so the bad pixels shifted out of alignment with one an- 
other. The result has no bad pixels and was calculated without 
interpolation. The mask sum gives the number of good pixels 
contributing to the average at each pixel. 

To apply the template, we shifted it according to Saturn's 
predicted position (see above) in each occultation frame, mul- 
tiplied it by Rhea's flux in that frame, down-sampled to the 
original spatial resolution, subtracted from the calibrated stel- 
lar frame, up-sampled the result to the scaled resolution using 
nearest-neighbor sampling, and performed aperture photome- 

^ http://ssd.jpl.nasa.gov/horizons.html 
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TABLE 2 
Observation and Light-curve Parameters 



Description 


Value 


Comment 


Telescope 


IRTF 




Instrument 


NSFCAM 




Image scale 


Of' 30 pixel- 1 




Wavelength, A 


2.28 ^m 


"Spencer 2.3" filter 


Bandpass, FWHM 


0.17 /im 




Time between frame starts 


0.2128 s 




Exposure time 


0.1276 s 


Dead time of 0.0852 s, see flie text 


Readout time 


0.009469 s 




Non-destructive reads (NDR) 


8 




Gain 


10 e" DN-1 




Read noise 


19 e" 


55 e- read' / VNDR 


Frame mid-times: 






First frame 


11:05:02.622 UTC 




First frame in upper baseline average 


11:16:44.939 UTC 




First frame in lower baseline average 


11:26:19.562 UTC 


No evidence of star in frames 


Last frame 


11:56:07.065 UTC 




Number of frames 


14,400 




Number of frames per baseline average 


300 




Photometry (upper baseline interval): 






Star full flux 


5439 DN 


Above lower baseline 


Rhea full flux 


5216 DN 


Above lower baseline 


Sky box 


12.55 DNpixer' 




Variability (upper baseline interval): 




Standard deviation 


Star 


2.2 % 


Scintillation uncorrected 


Star 


1.5% 


Scintillation corrected 


Sky box 


3.8% 




S/N per frame 


66 




S/N per scale height 


267 





try using the predicted stellar position. 

The 0.15-pixel centering uncertainty would lead to slight 
misalignments of the template. In the region of the limb, the 
template's gradient is typically 0.0018 DN per subpixel (the 
peak is 0.0054 DN per subpixel). Both these values are found 
after multiplying by Rhea's mean flux from Table 2 and di- 
viding by the square of the scaling (100, so that photometry 
yields an unsealed DN value like those in Table 2). The aper- 
ture boundary touches the limb gradient in only two places, 
however; so template shift errors contribute minimal noise. 
Shifting the template by one or two subpixels horizontally or 
vertically, subtracting it from the unshifted template, and per- 
forming aperture photometry yields -1.26 and -2.10 DN for 
the horizontal shifts and 1.21 and 2.08 DN for the vertical 
shifts, respectively. The photon noise from the unocculted 
star is 23 DN, from the sky in the aperture is 12 DN, and from 
other sources is 15 DN based on aperture photometry of the 
template. Template misalignment is thus a minor noise source 
even compared to photon noise. 

The template-subtracted stellar frames were quite flat ex- 
cept in the corner that included some ring light (see Figure 
1). Pixels in this region were in the sky annulus but not in the 
photometry aperture, and were efficiently rejected by taking 
the median of all good pixels in the sky annulus. Sky pho- 
tometry for Rhea also used the median, for consistency, as did 
that in the preceding paragraph. 

The stellar and Rhea photometry apertures both had di- 
ameters of 3 "6, four times the seeing FWHM. The abutting 
sky annuli had outer diameters of 7f'2 (see Figure 1). The 
photometry for the unocculted star and for Rhea closely fol- 
lowed each other, indicating that the main noise contributor 
was scintillation. We thus divided the stellar flux by Rhea's 
flux in each frame (just as we adjusted the template's flux). 



The light curve's upper and lower baseline levels are the 
medians of sections of the data before and after immersion 
(see Table 2). We used the median to reduce the effect of 
non-Gaussian residual scintillation. The lower baseline esti- 
mate and the template come from the same frames. Verify- 
ing that there is no star in the lower-baseline frames is cru- 
cial to the fits and inversions that follow. The lower baseline 
frames were aligned at the nominal Saturn position and added 
to make the template. Due to increasing ray bending, the nom- 
inal stellar position is a relatively stationary point on the limb 
at this time in the occultation. The template shows no evi- 
dence of the star and is smooth (see above). After template 
subtraction and scintillation correction, the median integrated 
value for the star aperture in the lower-baseline frames was 
(1.5 ± 6.0)x 10 '^ units of the full stellar flux over the level in 
the sky annulus. The mean was less than twice this value, and 
the uncertainty is the standard deviation of the mean. The star 
is effectively gone. 

Figure 2 presents the light curve. Its per-frame signal-to- 
noise ratio (S/N) is 66 for the unocculted star The French 
et al. (1978) atmospheric occultation noise parameter e^ = 
e^v^/H = 0.00375, where e is the standard deviation of the 
unocculted stellar flux in one-second bins divided by its mean, 
v^ is given in Table 1, and H, the scale height, is given in Ta- 
ble 1 of Paper 11 (baselines fixed case). The reciprocal of e^ 
is the S/N per scale height, and both are commonly used as 
figures of merit for occultation light curves. Our value of 267 
exceeds even that of the extremely bright 28 Sgr occultation. 
For example, the IRTF 28 Sgr immersion S/N per scale height 
was 192, using the fractional flux standard deviation per frame 
of 1.7% and exposure time of 0.25 s reported by Harrington 
et al. (1993) and ///v_l = 2.66 reported by Hubbard et al. 
(1997). The raw S/N for these occultations is limited by scin- 
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Fig. 2.— Light Curve for occultadon of GSC 0622-00345 by Saturn, observed at the NASA IRTF on 1998 November 14 UTC. The light curve has stable 
upper and lower baselines (dashed lines) and quite low photometric and terrestrial atmospheric scintillation noise. This makes it possible to determine the thermal 
profile of the atmosphere by numerical inversion. The light curve is characterized by many sharp, narrow spikes, due to small-scale refractive focusing by density 
fluctuations in Saturn's atmosphere. In some cases, these spikes have an amplitude of several times the ambient intensity in the light curve. The smooth, solid 
line is an isothermal model fit with baselines fixed (see Table 1 of Paper II for parameters), which we note in the text does not fit well near the lower baseline. 
Below the light curve, we mark the location of the half-light point (z = km) and the extrema of the reliable region of the inversion, using altitudes determined 
by the inversion. 



tillation rather than by photon statistics, eliminating 28 Sgr's 
8 mag brightness advantage over the event reported here. Our 
sHghtly slower event, improvements in infrared array technol- 
ogy, and the analysis methods described herein account for 
the rest. The electronic supplement to this article includes the 
light curve data file. 

A histogram of our upper baseline residuals closely follows 
a Gaussian distribution to ±3cr, but at >+3a it is slightly 
above the Gaussian, perhaps because scintillation produces 
occasional, transient flashes. The standard deviation in the 
difference between adjacent-frame fluxes is >93% of that ex- 
pected for uncorrelated errors, and reaches ^99% over an in- 
terval of 1 s. The residual correlation may come from the < 5 
ms interval between reads of the star and Rhea pixels. The 
power spectrum of the noise is thus still slightly stronger at 
lower wavenumbers than at higher wavenumbers (see below). 
In the occultation itself, one would expect correlated structure 
on the scale of the projected stellar diameter or the Fresnel 
diffraction scale (see Table 1), but both of these are below the 
2.7-km sky-plane resolution of the light curve. 

We computed the nominal occultation track from the 
ephemeris, and then offset by an additional 500 km E and 500 
km N, relative to Saturn's center, so that the absolute radius 
of the half-light level matched the oblate half-light surface 
determined by Hubbard et al. (1997) from 28 Sgr occultation 
observations. This shift is within the astrometric accuracy of 
the star relative to Saturn. 

Since the PSFs of Rhea and the star were nearly the same, 
we attempted to perform optimal stellar photometry, analo- 



gous to the optimal spectral extraction of Home (1986) and 
others. One of us (J.H.) used this method for exoplanetary 
secondary eclipses from the Spitzer Space Telescope (Deming 
et al. 2005), but developed it initially for the present analysis. 
We divided each background-subtracted point-source image 
by a normalized PSF template. Each pixel in the subarray thus 
became an estimate of the total flux of the point source, with a 
variance that increased (rapidly) with distance from the stellar 
centroid. We then computed an error-weighted (i.e., optimal) 
mean of all the good pixels in the frame. 

In principle, this process can reduce the noise by up to 70%. 
The optimal weighting discounts pixels far from the stellar 
centroid, and thus reduces the effect of read noise and residual 
photon and scintillation noise from scattered planetary light. 
There is also no photometry aperture size decision to make. 
However, the method requires an accurate PSF model and 
good centroids. The highly varying and non-Gaussian PSFs 
of these short-exposure images precluded any smoothing or 
analytic PSF model fit, so the PSF templates were just the 
background-subtracted, normalized images of Rhea, shifted 
to the centroided position of the star in its images and down- 
sampled to the original resolution. For both Rhea and the star, 
we used the background level from aperture photometry. 

Regrettably, while this procedure yielded a lower baseline 
that was three times less noisy than in the aperture-derived 
light curve, the upper baseline was about four times nois- 
ier. Close examination of the PSFs showed that noise and the 
fractional-pixel centroiding errors were sufficient to compro- 
mise this method. We decided, for the sake of consistency, not 
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Fig. 3. — Portion of the LED timing liglit curve and tlie best-fit model. 
Traces for botli the LED light curve and the model overplot each other at the 
top of the plot. The lowest trace presents the residuals, offset by -0.05 for 
clarity. 

to merge the initial aperture photometry with the final optimal 
photometry. However, a much brighter PSF standard might 
have substantially improved our photometry. Obtaining such 
a standard may be difficult in occultation experiments with 
bright stars, but we note that large satellites could still play 
this role in occultations of fainter stars, potentially increasing 
the number of useful occultation events. If optimal photom- 
etry were successful, the disappearance of scintillation noise 
with atmospheric depth would argue for the application of sig- 
nificance tests similar to the scaled noise tests of Paper 11, as 
opposed to the unsealed set. 

We did not acquire useful egress data. Atmospheric emer- 
sion occurred in the rings (see Figure 1), so, following Har- 
rington et al. (1993), we employed a 3.34 /im filter (wave- 
length measured warm) in which both planet and rings are 
dark due to overlapping absorption bands of water ice and 
methane. To increase the frame rate, NSFCAM's movieburst 
mode cannot display data during acquisition. It was thus not 
possible to see until afterward that the frames were empty. 
Possible reasons include a rise in the telluric water column 
(which would absorb all the light at this wavelength), a filter 
wheel problem, or (unlikely) a pointing problem. Two pro- 
cesses using shared memory on a dual-CPU camera computer 
could implement a movieburst mode with simultaneous image 
display without compromising the frame rate. 

3. TIMING AND ARRAY CLOCKING 

Relating each flux measurement to a location in the atmo- 
sphere requires precise knowledge of the time of each expo- 
sure. At the start of the first movieburst frame, NSFCAM's 
computer reads the time from a Global Positioning System 
(GPS) receiver The observer must find the delay between the 
GPS read and the actual first collection of photons, the actual 
exposure time, and the time to read the array (which is done 
many times per frame). For this purpose, there is a flashing 
light-emitting diode (LED) that is driven by the GPS and lo- 
cated in the camera. The LED turns on for 0.1 s at the start of 
each Coordinated Universal Time (UTC) second. It is often lit 
for only part of an exposure, resulting in count levels that are 
between the fully lit and unlit levels. By observing this LED 
with exactly the same settings as for the occultation, one can 
model the count levels to determine the timing parameters. 

We made such measurements before and after the occulta- 



tion. For each frame, we summed the flux in the sky subarray 
(which was lit by the LED) after subtracting from each pixel 
the average per-pixel flux in the star subarray (which was not). 

Our model starts with five lists of event times: frame re- 
sets, pedestal non-destructive reads (NDRs, see below), read- 
out NDRs, LED "on" times, and LED "off" times. It concate- 
nates the lists, sorts them in time, and steps through the re- 
sult to calculate the flux per frame. It thus takes into account 
the changes in effective sensitivity induced by each NDR. The 
five free parameters are delay after GPS read, integration time, 
readout time, and the two LED flux levels. We iteratively 
fit the model to the LED data by minimizing the rms residu- 
als. The relevant results appear in Figure 3 and Table 2. The 
electronic supplement to this article includes code that imple- 
ments the model in the Interactive Data Language (a product 
of ITT Visualization Systems, Boulder, CO, USA). 

The per-image timing uncertainty ("jitter") is difficult to es- 
timate, as fluctuations in the LED brightness, read noise, etc. 
all masquerade as jitter This parameter is important in stud- 
ies of detailed light curve structure, such as the shape and 
location of spikes, but this is not the focus of our present in- 
vestigation. 

Future instruments will likely read the GPS for each frame. 
However, an end-to-end test such as that provided by NS- 
FCAM's LED builds confidence that these increasingly- 
sophisticated instruments have consistent timing. It also 
quantifies the actual integration time and the delay between 
the GPS read and the first detected photons. 

We used NSFCAM's default "array reset clocking, double 
sample" mode. This mode performs a series of NDRs both 
before and after the exposure, and records the difference be- 
tween the two sets' averages. The multiple NDRs suppress 
read noise but reduce sensitivity for part of the exposure. We 
used the NSFCAM default of eight NDRs, but should have 
reduced dead time and increased read noise by using fewer 
NDRs. There were two effects of our long dead time. First, 
the photon noise was higher than it could have been. Since 
scintillation still dominated the noise even after our correc- 
tion, this is regrettable but not serious. Second, flux varia- 
tions occurring during either set of NDRs were only fraction- 
ally recorded, depending on when they occurred in the set. 
We may have missed some of the light in fast spikes. These 
correspond to the smallest-scale features in the inversion pre- 
sented in Paper II. Since the small scales have the smallest 
amplitudes (as expected, see Wavelet Analysis, Paper II), we 
conclude that our error did not strongly affect the results. 

4. CONCLUSIONS 

We have presented a light curve based on IRTF observations 
of the 1998 November 14 occultation of GSC 0622-00345 by 
Saturn. The light curve has a per-frame S/N of 66 and a per- 
scale-height S/N of 267, placing it among the best ground- 
based atmospheric light curves. We have described the obser- 
vational and analytic methods used to derive the light curve, 
in some cases for the first time in the literature, although 
these techniques are known among experienced occultation 
observers. Methods that improved the photometry include 
using Rhea as a scintillation and pointing standard, spatially 
subsampling the data by a factor of 10 for accurate aperture 
placement, and removing a template of scattered light. Sep- 
arate, identical, measurements of a flashing LED placed in 
the telescope beam provided an independent measure of the 
frame rate, when fit by a detailed model of the chip's sen- 
sitivity throughout the imaging cycle. Optimal photometry 
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shows promise at reducing the photometric noise level as the 
star dims, but requires a PSF standard that is substantially 
brighter than the occultation star. Occultation cameras with 
wider fields of view (but still able to operate with accurate 
timing at high frame rates) would facilitate finding bright PSF 
standards. 
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